Office  of  Naval  Research 


Q 

Contract  N50RI-76 


- - . - dS 

Task  Order  No.  1 ♦ NR -071 -012 


SOME  SPIN  WAVE  PROPERTIES  OF 
FERRIMAGNETIC  AND  ANTIFERROMAGNETIC 
SIMPLE  CUBIC  CRYSTALS 


BU 


J.  S.  Kouvel  and  Harvey  Brooks 


May  20, 1954 

Technical  Report  No.  198 

o 1 — - " — — " - ~~  — ' = ‘ *0 

Cruft  Laboratory 
Harvard  University 
Cambridge,  Massachusetts 


THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOD  DIRECTIVE  5200,20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELKA'JEj 
DISTRIBUTION  UNLIMITED, 


r 


Office  of  Naval  Research 
Contract  Nbori-76 
Task  Order  No.  1 
NR-071-01 1 


Technical  Report 
on 

Some  Spin  Wave  Properties  of  Ferrimagnetic  and 
Antifeirromagnetic  Simple  Cubic  Crystals 

by 

J.  S.  Kouvel  and  Harvey  Brooks 


May  20,  1954 


The  research  reported  in  this  document  was  made  possible 
through  support  extended  Cruft  Laboratory,  Harvard  University, 
jointly  by  the  Navy  Department  (Office  of  Naval  Research),  the 
Signal  Corps  of  the  U.  S.  Army,  and  the  U.  S.  Air  Force,  under 
ONR  Contract  N5ori-76,  T.  O.  1. 


i 

Technical  Report  No.  198 


v 


Cruft  Laboratory 
Harvard  University 
Cambridge,  Massachusetts 


C 


r 


TR  198 


Some  Spin  Wave.  Properties  of  Fer rimagnetic  and 
Antiferromagnetic  Simple  Cubic  Crystals 

by 

J.  S.  Kouvel  and  Harvey  Brooks 
Cruft  Laboratory,  Harvard  University 
Cambridge,  Massachusetts 


1 . Introductory  Remarks 

Several  theories  on  the  electronic  spin  resonance^  in  ferrimagnetic 
and  antiferroma  gnetic  media  have  been  recently  developed  to  a point  where 
their  results  may  be  readily  compared  with  experiment.  These  have  all 
been  based  on  simple  semiclas sical  principles,  and  hence,  must  await  con- 
firmation from  more  rigorous,  quantum  mechanical  work,  also  now  in 
progress.  ^ The  main  purpose  of  this  paper  is  a semiclassical  spin  wave 
treatment  of  certain  "bulk"  properties  (namely,  the  specific  heat  and  satura- 
tion magnetization  at  low  temperatures)  of  a simple  cubic  ferrimagnetic 
or  antiferromagnetic  crystal;  vhe  results,  we  hope,  will  also  be  amenable 
to  at  least  qualitative  comparison  with  experiment.  For  the  sake  of  sim- 
plicity, this  study  is  confined  to  a simple  cubic  single  crystal  which  we 
postulate  to  have  a single  axis  of  magnetic  anisotropy  (along  a cube  edge) 
instead  of  the  cubic  anisotropy  generally  associated  with  such  a crystal. 

The  effects  of  a uniform  magnetic  field  applied  along  the  anisotropy  axis 
are  studied  lor  both  ferrimagnetic  and  antiferromagnetic  cases;  moreover, 
for  the  latter  case,  the  situation  in  which  the  field  is  applied  perpendicular 
to  the  anisotropy  axis  is  also  considered. 

Since  our  results  can  be  readily  made  to  apply  to  the  ferromagnetic 
case  of  the  simple  cubic  model,  v/e  are  able  to  check  them  directly  with 
the  results  of  previous  spin  wave  analyses  of  the  ferromagnetic  problem . 
Consequently,  it  becomes  possible  to  trace  the  similarities  and  differences 
of  the  "bulk"  spin  wave  properties  of  all  three  cases:  ferromagnetic,  ferri- 
magnetic, and  antiferromagnetic. 


c 
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2.  The  Ferrimagnetic  Case 

The  atomic  arrangement  in  our  simple  cubic  model  is  ordered  with 

respect  to  the  net  electronic  spin  per  atom  in  the  manner  shown  in  Figure 

1.  The  ground  state  of  this  system  is  assumed  to  be  such  that  the  electronic 

spin  of  each  atom  of  sublattice  (1)  is  directed  upwards  (in  the  z direction) 

while  the  spin  of  each  atom  of  sublattice  (2)  is  directed  downwards.  Though 

there  is  no  rigorous  theoretical  justification  for  this  assumption,  it  is  hoped 

that  the  arguments  of  Anderson  and  others2  for  very  nearly  this  type  of 

long-range  order  in  the  ground  state  of  a simple  cubic  antiferromagnetic 

may  be  also  applicable  in  the  fer rimagnetic  case.  In  fact,  these  arguments, 

being  basically  classical,  should  have  increased  validity  for  large  atomic 

spins,  which  is  usually  the  case  for  ferrimagnetics . Furthermore,  Ne'el's 

3 

phenomenological  theory  of  the  saturation  magnetization  of  ferrites, which 
has  enjoyed  remarkable  agreement  with  experiment,  is  firmly  baaed  on  an 
antiparallel  alignment  of  unequal  atomic  spins. 

The  energy  of  negative  exchange  interaction  which  is  responsible  for 
this  antiparallel  alignment  may  be  expressed  as 


2 J ) S.  • S . 

L—s  -i  -j 

i.j 

so  that  J,  the  negative  of  the  exchange  integral,  is  a positive  quantity.  For 
the  present,  we  shall  consider  only  nearest  neighbor  exchange  interactions; 
the  additional  effects  of  next-to-nearest  neighbor  interactions  are  discussed 
in  Appendix  A. 

An  axial  magnetic  anisotropy  which  tends  to  direct  the  spins  along  the 
z-axis,  is  assumed  to  have  an  energy  associated  with  it  as  follows: 


EK  .S.  sin2  8 . 
J J J 


j 


or 


(S  ,+S  .) 

-*j  -yj 


where  K.  is  the  anisotropy  energy  per  unit  spin  of  S., 
J J 

S.  makes  with  the  z-axis. 

-J 


and  6 . is  the  an 
J 


gle 


Furthermore,  the  syst  - i is  subjected  to  a uniform  magnetic  field, 
H,  in  the  z direction,  which  ^ives  rise  to  the  energy, 
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-H 


S. 

-J 


where  (3  is  the  Bohr  magneton  and  gj  is  the  spectroscopic  splitting  factor 
for  the  electronic  magnetic  moment  of  the  jth  atom. 

The  sum  of  these  three  energy  terms  is  taken  to  be  the  total HamiLtonian 

of  the  system.  Throughout  this  report,  we  have  thus  neglected  the  effects 

of  magnetic  dipole -dipole  interactions  between  spins,  except  insofar  as  they 

4 

are  included  phenomenologically  in  the  anisotropy  energy;  for  a structure 
having  the  simple  cubic  symmetry  shown  in  Figure  1,  they  may  be  expected 
to  be  small. 

The  equation  of  motion  of  tne  4th  spin  vector  may  be  written  as 


2K 


t 


ST-  5-1  x (-x£+5y£'  T * EL 


) + gp  PS  p X H 


(2.  1) 


Indicating  the  two  sublattices  by  superscripts  (1)  and  (2),  we  express  the 
spin  vectors  in  their  rectangular  components  as  follows: 


:(D 


A 

= 3 


S^+j  S 

XX  Jy  y 


w+j 

\r  J 


si 

z 1 


(2.  2a) 


4 S(2)4  S(2)4  S 
■^x  x Jy^y  ^z^2 


(2. 2b) 


The  z components  of  spin,  to  a first  approximation,  are  considered  to  be 
invariant  in  time  and  space.  This  is  equivalent  to  neglecting  terms  quad- 
ratic in  the  S’s  and  S ’s.  Furthermore,  we  assume  that  all  x and  y spin 
x y 

components  are  expressible  in  the  form 


;(n) 


:(n) 


S'  '(x.y.z.t)  = S'  (q  ,a  ,a  ,t)  exp 


^2-rri  (a. 


x + a y + a z)/a 
x y7  z 


(2.3) 


where  a is  the  distance  between  nearest  neighbor  atoms  and  a +a  +a  - a 

° —x  — y — z — 

is  a dimensionless  phase  vector  (corresponding  to  Anderson's2  X./2tt  and 

5 ~ 

to  Herring  and  Kittel’s  ka/2Ti).  The  equations  of  motion  for  the  x and  y 

spin  components  are  then  found  to  be 


r 
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ks(1)  = 1 2JS  . y S(2)+  ( 12JS-  + 2I<  + g.'PHlS*1’  (2.4a) 

x j.  q>  y dd  a a y 

W2)= -12JS2^isJ1)+(-12JS1  - 2K2+  g2pH)S*2)  (2.4b) 

ks^1^ -12JS1vaS^2)+(-12JS2-  2K:  - g^HjS^  (2.4c) 

ks(2)  = 1 2JS , y S(1)+  ( 12JS  . + 2K,  - g,pH)si,2)  (2.4d) 

y Ld  CL  X A dd  Ld  x 


where 

y = (cos  2im  + cos  2iTa  + cos  2ira  )/3 
*a  x y z 


Wc  now  assume  that  the  time  dependence  of  the  x and  y spin  components 
is  of  the  form 


By  means 
yield 


i<n)+ 

X 

iS(n)=  s^e-*** 

y + 

(2.5a) 

in)- 

iS^n)=  sin)eiwt  . 

(2.5b) 

of  (2.5a),  (2.4a)  is  combined  with  (2.4c),  and  (2.4b)  with  (2.4d)  to 


kosj1^  12JSrYaS|2)+  (12JS2  + 2*^  + g^Hjsj1*  (2.6a) 

kws|2)= -12JS2vqs|1)+(-12JS1  - 2K2+  g2pH)sj2)  (2.6b) 


[(2.5b)  leads  to  the  same  equations  with  the  subscripts  + replaced  by  -]  from 
which  the  following  solutions  for  the  frequencies  are  immediately  obtained: 

kw  = -6J(SrS2)  + Kl  - K2  + (gl+g2)pH/2 
±[(6J(SrSj}2 

+ [K1+K2+(g1-g2)pH/2j  [l2J(S1+S2)+K1+K2r(g1-g2)pH/2] 

2 2 ~\  1/2 

+ (12Jr  S^l-y^)  (2.7) 

At  very  low  temperatures,  where  the  spin  wave  theory  is  most  likely 
to  be  valid,  only  the  spin  waves  of  very  long  wavelength  compared  to  the 


FIGURE  I 


t 


TR  198 


-5- 


lattice  spacing  may  be  expected  to  be  excited  appreciably.  Hence,  a ' /ill 
range  over  very  small  values  and 


Y 


o/  1 - 

rt 


or 


(2.8) 


Moreover,  for  a ferrimagnetic  with  a reasonably  high  Curie  temperature, 
bJtS^-S^)  is  substantially  larger  than  the  terms  in  (2.7)  that  contain  the 
anisotropy  constants,  the  applied  field,  or  a.  It  follows  that  the  frequencies 
of  the  two  modes  are  very  different  from  each  other.  We  can  ignore  the 
higher  or  "optical"  mode  since  its  contributions  to  the  specific  heat  and  to 
the  temperature  change  of  magnetization  are  negligible  at  temperatures 
well  below  the  Curie  point,  and  consider  only  the  lower  or  "acoustical"  mode 
(corresponding  to  the  choice  of  the  positive  sign  in  (2.7)).  It  should  be 
noted  that  a discrepancy  between  the  g-factors  of  the  two  sublattices  causes 
an  apparent  field-dependent  increase  of  the  anisotropy;  this  was  recently 
discussed  by  Wangsness^  in  his  analysis  of  ferrimagnetic  resonance.  Con- 
versely, a difference  between  the  two  anisotropy  constants  shows  up  as  an 
apparent  increase  of  the  g-factors  (or  of  the  external  field).  Hence,  keep- 
ing the  possibility  of  these  effects  in  mind,  we  proceed  to  simplify  (2.7)  by 
setting  = K and  g^  = g^  = g.  We  then  expand  the  square  root  in 

(2.7)  and  obtain  for  the  "acoustical"  mode 


Itw  di  gpH+2K 


fS1+S.+K/'6Jl 


r 


12JSlS2^Z  2^ 

-spsr  \T  a ) ' 


(2.9) 


For  comparison,  the  spin  vibrational  frequencies  for  a ferromagnetic 
simple  cubic  structure  may  be  readily  obtained  from  (2.7)  by  setting 
Sj  = -S^  = S,  Kj  = -K^  = K,gj  = g ^ = §>  and  “J  = J'»  in  which  case  it  is 
founa  that 

h<o -gpH  + 2K  + 6J'S  a2^  . (2.10) 

Thus,  except  for  the  generally  small  effects  of  anisotropy  and  applied 
magnetic  fields , cj  is  proportional  to  a in  both  the  ferrimagnetic  and 
ferromagnetic  cases.  Such  a relationship  between  w and  a,  as  has  been 
shown  in  previous  discussions  of  the  ferromagnetic  case,  leads  directly 


I 

w 


r 
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to  a T law  of  the  specific  heat  and  of  the  decrease  of  saturation  magnetiza- 
tion. We  shall  now  make  use  of  equation  (2.9)  to  derive  the  actual  expressions 
for  these  properties  of  a ferrimagnetic . 

Quantizing  the  normal  modes  of  our  classical  spin  wave  problem,  we  are 

7 

led  to  the  following  expression  for  the  specific  heat  per  atom: 


°o  , H - 

r - 1-  C ->  e (4™ 

C a ■ k / TTT 

•4=o  <e  -1' 


da) 


(2.11) 


where  5 = ^iw/kT . By  means  of  (2.9),  (2.  11)  is  transformed  to 


C.  = 

a 


4ir 


(SrSz)kT 


-^57 


3/2  xo 


J 


5H' 


(e5-  1)Z 


1/2 


dB 


(2. 12) 


in  which  = gpH'/kT  where  H' , the  total  equivalent  magnetic  field  exerted 
on  the  3yatem,  is  defined  a« 


H'  = H + 


2K  (VS2+ K/6JN 


gp  T^5-2“ 


(2.13) 

co 


■c  'p  _2  y— t _nt 

In  evaluating  the  integral  in  (2.12),  we  first  express  e3  (e  -1)  as  / ,ne 


and  then  substitute  u = (5-5^,)  ^2.  The  integralthus  becomes 


n=  1 


00  w 00 

XI  2re  J u2(u4  + 25h,u2  + 5h^) 

n=  1 o 


-nu 


8 

From  tables,  we  find  that 


du 


00 


,a+ 1 a xp 

2 p 


(2. 14) 


and  the  integral  is  thus  finally  determined  as 

00 


15\|rr 

8 


XI  +Fn§H'+TT(n^H')2jn 


-5/2  _n^H' 
e 


For  the  special  case  of  = 0,  this  expression  reduces  to  simply 
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15  Vir 
8 


or 


0/2) 


where  the  Riemann  zeta  function,  £(5/2)  = 1.341.  Consequently,  the  spin 
wave  contribution  to  the  specific  heat  of  the  ferrimagnetic  system  may  be 
written  as 


(Ca^H--° 

(ca^H;=o 


1 + 5 n5jp  + it; 


15 


3H' 


-5/2 
n e 


-n^ 


H' 


(2.  15a) 


and 


<Ca>5 H.=» 


= 0.113k 


(S.-SJkT' 


3/2 


(2.  15b) 


if  the  anisotropy  and  applied  magnetic  fields  are  negligibly  small,  the 

3/2 

specific  heat  follows  a T law  exactly.  However,  should  they  be  ap- 
preciable, (2.  15a)  indicates  that  the  temperature  dependence  of  the  spe- 
cific heat  no  longer  follows  such  a simple  relationship.  In  fact,  from 
Figure  2a,  where  the  right  hand  side  of  (2.  15a)  is  shown  as  a function  of 
kT/gpH'  (or  1/5jji )*  it  Is  evident  that  at  sufficiently  low  temperatures,  any 
finite  anisotropy  or  applied  field  will  ultimately  cause  the  temperature 
dependence  of  the  specific  heat  to  be  exponential.  To  illustrate  the  change 
of  the  temperature  dependence  of  the  specific  heat  with  applied  field,  we 
have  plotted  in  Figure  2b  the  specific  heat  against  absolute  temperature 
f°r  g pH'/k=0,  1 , and  2 (i . e . , for  H'=  0,  7430,  and  14860  gauss , if  it  as  - 
sumed  that  g = 2) . Finally,  in  Figure  3,  we  have  plotted  the  right  hand 
side  of  equation  (2.  15a)  against  tc  show  how  the  specific  heat  varies 
with  applied  field  at  constant  temperature.  The  monotonic  decrease  of 
the  specific  heat  with  increasing  H may  be  understood  qualitatively  from 
the  fact  that  the  spin  excitational  energies,  liu,  increase  with  H (see  equa- 
tion (2.9)),  thus  contributing  less  to  the  specific  heat. 


If  the  zero-point  energy  is  neglected  for  reasons  analogous  to  those 
2 

given  by  Anderson  for  an  antiferromagnetic  system  of  large  atomic  spin, 
the  total  ener  gy  of  our  ferrimagnetic  system  may  be  written  as 


r 
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W = -6NJS  - NgpHtSj-S^/Z  + ^ nQkw 

a 

3 

where  n^,  a function  of  a,  is  the  integral  occupation  number , and  N = 1/a  , 
the  number  of  atoms  per  unit  volume  of  the  system.  It  follows  that  the 
partition  function  is 

Z = exp  | ( 12N JS  + NgPH(S,-SE))/2kT]  ] [ ( 1 - e"^)"1  . 

i Q 

Using  this  expression  for  Z together  with  (2.9),  we  determine  the  satura- 
tion magnetization  as 

M = kT  g|j(lnZ)  = Ngp(SrS2)/2  - g p (e5  - 1)  ‘ 1 . 

a 

Thus,  Mq>  the  saturation  magnetization  at  absolute  zero,  is  Ngj3(S  ^ -S2)/2, 
and  the  variation  of  the  relative  saturation  magnetization  is  found  to  be 


M -M  AM 
o o 

— n or  n — 


o 2 \ ’ , ^ , *-l  2 C 4ira  da  i*  it \ 

Njspsp  I,1'  - 11  Vs!  J „~X7  • ' ■ ' 


By  means  of  (2.9),  we  convert  (2.  16)  to 


A M 

o 2ir  (S 


1 \ (Sj-SjIkTj  3/2  r”  (S-5H,)1/2dS 


(2.17) 


t _ ^ oo  t 

In  determining  the  integral  in  (2.  17),  we  first  express(e  - 1)  as  ^ e n~and 
use  the  substitution,  u = Then,  by  the  application  n 

of  (2.  14),  we  get  for  this  integral 


(v?/2)  n_3/z 


which,  ior^^,  = 0,  reduces  to 
n. 


u 


FIGURE  2b 


A,CaAND  A2A0M/M0 
FOR  H1  =0 
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(V5r/2)  £(  3/2)  or  (Vir/2)  (2.612)  . 


Hence,  the  relative  saturation  magnetization  of  the  ferrimagnetic  structure 
may  be  written  as 


(2. 18a) 


(2.  18b) 


3/2 

Similar  to  the  specific  heat,  A Ivl/M  departs  from  a T~  law  when  the 

o o 

anisotropy  and/or  applied  magnetic  fields  are  sufficiently  large.  The  ap- 
proach to  an  essentially  exponential  temperature  dependence  of  AqM/Mo  at 
low  temperatures  is  shown  in  Figure  2a  where  the  right  hand  side  of  equa- 
tion (2.  18a)  has  been  plotted  against  kT/gpH' . In  Figure  2b,  AqM/Mo  has 
been  plotted  against  absolute  temperature  for  gfiH'/k  = 0,  1,  and  2.  It  is 
quite  evident  that  AqM/Mq  is  much  more  strongly  influenced  than  the 
specific  heat  by  anisotropy  and  applied  fields.  This  is  also  obvious  from 
Figure  3 where  the  right  hand  side  of  equation  (2.  18a)  has  been  plotted  as 
a function  of^jj,. 

We  previously  showed  that  the  expressions  for  the  spin  excitational 
energies  for  both  the  ferrimagnetic  and  ferromagnetic  cases  of  the  simple 
cubic  model  have  the  same  form  (see  equations  (2.9)  and  (2.  ]0)).  If  we 
now  go  through  the  same  procedure  for  the  ferromagnetic  case  by  start- 
ing with  (2.  10)  instead  of  (2.9)  and  letting 


H'  = H + 2K/g  (3 


(2.19) 


we  find  that  equations  (2.  1 5a)  and  (2.  18a)  apply  exactly,  while  (2.  15b)  and 
(2.  18b)  should  be  replaced  by 


(Ca'SH,  = 0 


= 0.  113k  (kT/2J'S) 


3/2 


(2.20) 


and 


r 
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( A M/M  )w  _n=  (0.0587/S)  (kT/2J'S)3/2 
o o -5j^i  ~ u 


(2.21) 


9 

respectively.  Equation  (2.21)  is  in  perfect  agreement  with  Miller's  ex- 
1 0 3/2 

tension  of  Bloch's  original  T law,  for  the  general  case  of  S>  1/2;  Mott 
and  Jones  ^ have  obtained  (2.20)  for  the  case  of  S = 1/2. 


Holstein  and  Primakoff  have  investigated  quantum  mechanically 
the  field  dependence  of  the  intrinsic  magnetization  of  a ferromagnetic.  If 
their  magnetic  dipole -dipole  interaction  terms  are  neglected,  their  expres- 
sion that  is  analogous  to  equation  (2.  18a)  but  is  valid  only  for  small  values 
of  5jji  i may  be  written  as 


(2.22) 


and  is  represented  in  Figure  3 by  the  dashed  curve.  The  agreement  at 
low  values  of  5t»i  with  the  curve  corresponding  to  equation  (2.  18a)  is 
very  good.  More  important,  perhaps,  is  the  qualitative  fact  that  both  e qua 
tions  (2.  18a)  and  (2.22)  predict  a sharp  decrease  of  AqM/Mo  as  is  in- 
creased from  zero.  It  should  be  remembered,  however,  that^jji  cannot 
be  zero  if  the  anisotropy  of  either  the  ferrimagnetic  or  ferromagnetic 
structure  is  not  zero  (see  equations  (2.  13)  and  (2.  19)). 


3 . The  Antiferromagnetic  Case 

When  Sj  = = S in  the  spin  structure  shown  in  Figure  1 and  the 

g-factors  of  the  two  sublattices  are  assumed  equal,  the  net  saturation 
magnetization  of  the  system  is  zero.  Hence,  the  system  is  antiferro- 
magnetic. If  we  assume  further  that  the  anisotropy  constants  associated 
with  the  two  sublattices  are  equal,  the  spin  vibrational  frequencies  for 
this  antiterromagnetic  case  are  found  directly  from  (2.7)  to  be 

U = gpH  + [4K(12JS  + K)  + (12JS)2(4ir2a2/3)]  1/2  (3.  1) 

when  we  also  make  use  of  (2.8).  Consequently,  in  the  absence  of  anisot- 
ropy and  applied  magnetic  fields , w is  directly  proportional  to  a.  This 
particular  dispersion  law,  as  has  already  been  pointed  out  by  others,^ 

i 

/ 


I 


r 
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3 3/2 

gives  rise  to  a T variation  of  the  specific  heat  (as  compared  to  the  T 

variation  in  the  fer rirnagnetic  ferromagnetic  cases). 

The  transition  from  (2.9)  to  (3.  1)  for  the  frequencies  of  spin  vibra- 
tion is  a discontinuous  one  when  is  reduc  ed  from  1/2  to  C,  if,  as  may 

be  expected,  and  must  have  half-integral  values  in  order  that  a spin 
wave  analysis  be  valid.  Hence,  although  the  antiferromagnetic  case  rep- 
resents an  extreme  case  (i.e.  , when  S^S^)  for  the  ferrimagnetic  spin 
structure  of  Figure  1,  far  less  similarity  may  be  expected  between  the 
spin  wave  properties  of  an  antiferromagnetic  and  a ferrimagnetic  than  be- 
tween those  of  a ferrimagnetic  and  a ferromagnetic. 

The  existence  of  the  two  "acoustical”  modes,  corresponding  to  + and 

13 

- in  (3.1),  was  first  pointed  out  by  Kittel  in  connection  with  antiferro- 
magnetic resonance.  Equation  (3.  1),  itself,  was  derived  more  recently  by 
Keffer  et  al . ^ 

It  will  be  shown  later  that  there  is  a change  of  the  spin  state  from 

that  pictured  in  Figure  1 when  g(3H  exceeds  2[K(12JS-K)]  ^2,  which  is 

1/2 

smaller  than  2[K(12JS  + K)]  . Thus,  the  square  root  in  (3.  1),  for  the 

range  of  applied  fields  for  which  this  equation  is  valid,  will  always  be' 
larger  than  gpH,  and  we  should  keep  in  mind  that  it  is  the  absolute  magni- 
tude of  oi  which  must  be  inserted  in  (2.  11)  to  determine  the  specific  heat. 
Let  us  write  in  place  of  (3.  1), 

k|w|  = [4K(  1 2JS  + K)  + (12JS)Z  (4tt2qZ/3)]  1/2  + gpH  (3.2) 

Inserting  (3.2)  in  (2.  11),  we  obtain  the  following  equation  for  the  spinwave 
specific  heat  per  atom: 


.3/2  3 r"  1/2 

ca  = 'tbs’  J •*? 

GO 

G(£j)  = 52e^(e^  - 1)  2 = [('§'/2)/sinli  (^/2)]  2 =^^n5Z  e (3-3) 

n=  1 

where  B - nl  w 1/kT  , Bu  = 8 pH/kT , and  § = 2[K(  12JS  + K)1  ' /k'l  . For  negligible 

n o 


t 
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anisotropy  and  external  fields  (i . c . , for^Q  and  0),  the  integral  in 

(3.3)  is  readily  evaluated  by  the  use  of  the  summation  form  of  G(^)  and 

15 

the  tabulated  integral 
to 


-ax  , , . , b+  1 

e dx  = b ; /a 


and  the  final  result  for  the  specific  heat  is 


(C  ) 


aW° 


_ (3) 3/2  4!  t(4)k 

2tt2 


kT 

12JS 


f) 


3 

times 


2 


= 13.7k  (kT/1 2JS) 3 . (3.4) 

The  multiplication  by  2 is  required  to  account  for  the  two  modes  of  spin 

3 

vibration,  which  in  this  special  case  are  degenerate.  Thus,  a simple  T 

law  is  followed.  However,  in  the  more  general  case  of  finite  anisotropy 

(with  or  without  an  externally'  applied  field),  it  is  obvious  that  the  specific 

3 

heat  deviates  from  a T law  by  an  amount  determined  by  the  integral  in 
(3.  3).  Since  we  have  been  unable  to  evaluate  this  integral  in  its  general 
form  analytically  in  any  simple  manner,  we  have  resorted  to  graphical 
integration. 


Let  us  first  examine  the  situation  when  there  is  a finite  anisotropy 
but  no  applied  magnetic  field.  In  Figure  4a,  the  integrand  of  the  integral 
in  (3.  3)  has  been  plotted  against  5 for  various  values  of  ^ . These  curves 
were  then  integrated  graphically;  the  results , multiplied  by  the  coefficient 
shown  in  (3.  3)  (and  by  2 to  allow  for  the  two  degenerate  modes  of  spin 
vibration),  lead  directly  to  the  specific  heat  values  shown  plotted  against 
in  Figure  4b.  It  is  evident  from  this  drawing  that  a finite  magnetic 
anisotropy,  especially  since  ^ q depends  roughly'  on  the  geometric  mean 
of  the  anisotropy  and  exchange  energies  (the  latter  generally  being  by  far 
the  larger),  would  always  cause  a decrease  of  the  specific  heat.  The  lower 
the  temperature,  the  more  pronounced  will  be  this  effect.  From  Figure  5a 
where  the  specific  heat  normalized  by  its  no-anisotropy  value  is  shown  as 
a function  of  l/gQ,  we  should  note  that  at  sufficiently  low  temperatures  a 
finite  anisotropy  would  result  in  an  essentially  exponential  temperature 
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dependence  of  the.  specific  heat.  Figure  5b  may  serve  to  illustrate  more 

directly  the  effects  of  a finite  magnetic  anisotropy  on  the  specific  heat 

3 

and  its  temperature  dependence  (i.e.,  the  deviations  from  a pure  T^law) 

With  the  application  of  a magnetic  field  on  the  antiferromagnetic 
system,  the  degeneracy  in  the  spin  vibrational  energies  is  removed  (see 
equation  (3.2)),  and  each  of  the  two  modes  may  be  expected  to  contribute 
differently  to  the  specific  heat.  As  an  example,  we  have  set  = 1 and 
have  plotted  in  Figure  ba  the  integrand  in  equation  (3.  3)  against  *5  for 
several  values  of  ^ . For  § i 0,  the  two  parts  of  the  integrand,  corre- 

sponding  to  the  lower  energy  (L..F.)  mode  and  the  higher  energy  (H.E.) 
mode,  are  both  represented.  Integrating  these  curves  graphically,  v/e 
obtain  values  for  the  specific  heat  which,  in  Figure  feb,  have  been  plotted 
against  ^ It  is  clear  from  this  figure  +hat  for  an  increase  of  the  applied 
magnetic  field,  the  increase  of  the  contribution  of  the  Lower  energy  mode 
to  the  specific  heat  more  than  compensates  for  the  decrease  of  the  higher 
energy  mode  contribution.  The  net  result  is  an  increase  of  the  spin  wave 
specific  heat  of  the  antiferromagnetic  system.  This  contrasts  strongly 
with  the  fer rimagnetic  and  ferromagnetic  cases  for  which  the  spin  wave 
specific  heat  decreases  monotonically  for  an  increasing  applied  field  (see 
Figure  3) . 

By  a completely  classical  analysis  of  the  ground  state  of  the  atomic 
spins  of  a simple  cubic  antiferromagnetic  structure  (discussed  more  fully 
in  Appendix  B),  we  have  found  that  when  the  magnetic  field  applied  along 
the  anisotropy  axis  is  raised  above  a certain  critical  value,  an  abrupt 
transition  of  the  ground  state  spin  configuration  may  be  expected  to  take 
place.  In  Figure  7,  we  indicate  schematically  the  arrangement  of  the 
atomic  spin  vectors  of  the  two  sublattices  for  State  I,  the  original  state  of 
the  system  which  we  have  been  discussing,  and  for  State  II,  the  state  when 
the  applied  field  is  higher  than  the  critical  value  . It  may  be  noted  that  in 
State  II  the  spin  vectors  are  not  aligned  antiparallel  to  each  other,  and 
hence  give  rise  to  a net  magnetization.  In  addition,  at  the  point  of  transi- 
tion between  states,  there  is  a correspondingly  abrupt  change  of  themag- 

/ 1 h 

netic  susceptibility  parallel  to  the  applied  field,  as  Neel  ' was  th<»  first 
to  point  out . 
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Our  primary  interest  here,  however,  in  this  tiauSl  tion  of  the  spin 
state,  is  in  its  effects  on  the  specific  heat  of  the  system.  The  spin  wave 
problem  for  the  spin  configuration  of  State  II  must  therefore  be  solved. 
As  shown  in  Appendix  B,  tbs  condition  for  which  the  transition  between 
States  I and  II  occurs,  i« 

gpH  = 2[K(12JS  - K)]  1/2  , (3.5) 

and  for  State  II,  the  angle  0,  (see  Figure  7)  satisfies 

sin  0 = gpH/2{12JS  - K)  . (3.6) 

Thus,  if  the  field  applied  parallel  to  the  anisotropy  axis  is  increased 
beyond  its  value  in  (3.  5),  9 will  increase  monotonicallv  until  gf}H  = 
2(12JS  - K),  when  the  spins  of  both  sublattices  will  have  been  brought 
into  parallel  alignment  with  each  other  and  with  H.  However,  since  the 
magnitude  of  the  exchange  integral,  J,  of  any  antiferromagnetic  crystal 
whose  Curie  point  is  above  4°K,  is  probably  of  the  order  of  (10)  *^ergs 
or  more  (see  next  section),  it  would  not  be  possible  in  this  case  to 
achieve  an  appreciably  large  value  of  0 with  any  physically  realizable 
magnetic  field.  Nevertheless,  since  it  is  fairly  reasonable  to  expect 
that  K <f^l2JS  for  most  antiferromagnetic  crystals,  it  is  still  possible 
to  achieve  the  applied  fields  required  by  (3.5)  and  even  somewhat  higher 
fields . 


For  the  spin  configuration  of  State  II  as  shr-’n  in  Figure  7,  the 
equation  of  motion  of  a particular  spin  vector  m be  expressed  by 
equation  ( 2.  1),  but  the  rectangular  components  of  the  spin  vectors  must 
now  be  written  as 


i(1) 

S<2) 


= jx(S  cos  0 ^ S 

= j ( -S  cos  0 + S 
Jx 


A ( 1) 

+ j s'1' 

y y 


(2) 

X 


)+j  S(2) 

' Jy  y 


(S  sin  0 + S^l\ 
z 9 

(S  sin  9 + S^2)) 


(3.7a) 

(3.7b) 


the  dynamic  components  having  their  sublattice  designated  by  super- 
scripts. It  is  then  assumed  that  the  dynamic  components  have  the  spatial- 
wave  form  of  (2.3).  The  resulting  equations  of  motion  for  these  com- 
ponents are: 


z z 
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12JSyq  sin  6 S^2)  - [(12JS  - 2K)  sin  0 - gpHjS^  (3.8a) 

R2)=  12JSyq  sin  6S^l)  - [( 12JS  - 2K)  sin  0 - gpH]S^2;  (3.8b) 

nS^1^ -12JSYQSin0S(x2)  + [( 1 2 JS  - 2K)  sin  0 - gpHjS^ 

+ 12JSYqcos  0S^2)  + ( 1 2JS  - 2K)  cos  6 (3.8c) 

no(2)=  -IZJSy  sin  0S(1)  + [( 1 2JS  - 2K)  sin  0 - gpH]S(2) 
y ‘a  x L J x 

- lZJSYaCos  0 - ( 1 2JS  - 2K)  cos  6 S^2)  (3.8d) 

1s(1)=-12jSy  cos8S(2)  - 12JS  cos  8S(1)  (3.8e) 

z 'a  y y 

W2)  = 1 2JS y cos  0S(1)  + 12JS  cos  0S'2'  ( 3 . 8f ) 

z 'a  y y 


when  terms  higher  than  the  first  order  in  the  dynamic  components  are 
discarded.  Since  the  spin  vectors  may  be  thought  to  oscillate  about  their 
static  equilibrium  positions  defined  by  (3.6),  no  static  (zero  order)  terms 
appear  in  *he  above  equations.  From  a consideration  of  the  time  invariance 
of  the  spin  vector  magnitudes,  it  follows  from  (3.7)  that  the  x and  z dynamic 
spin  components  are  interrelated  by  the  equations, 

M)  cd)  c(2)  s(2) 

_ X =_^z  _q(l)  X _ z _ (2) 

sin  0 cos  0 ' u a sin  0 cos  0~  u ’ 

which  also  serve  to  define  a new,  convenient  set  of  spin  components,  ^ 

( 2)  u 

and  . Removing  0 from  equations  (3.8)  by  means  of  (3.6),  we  find  that 

in  terms  of  the  u and  y spin  components,  equations  (3.8a)  and  (3.8b)  are 
identical  with  (3.3e)  and  (3.8f),  respectively.  The  resultant  four  inde- 
pendent equations  are 


W^=  - 1 2JSy  S(2)  - 12JSS(1' 
u ’a  y y 


12JSy_  s 


(1) 


a y 


(3.9a) 


+ 12JSS(2) 

y 


(3.9b) 


r 
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kS(1)=  12JSY  (1  -R2/2)S(2) 
y 'a  H u 

+ [ 12JS  - 2K{1  - R^/4)] 


(3.9c) 


= -12JSy  (1  - R2 /2)  S( 
y V H u 


(0 


- [ 12JS  - 2K{1  - R2  /4)]  S(2' 
n u 


( 3. 9d) 


where  R^  = gpH/(12JS  - K).  Assuming  each  dynamic  spin  component  varies 

harmonically  with  time  (i.e.,  as  e~lwt),  we  proceed  to  combine  (3.9a)  with 

(3.9b),  and  (3.9c)  with  (3.9d)  to  form  two  pairs  of  linear  equations,  one  in 

S(1)-S(2)  a.-.d  S(1)+  S(2),  and  the  other  in  S(1)+S(2)  and  S(1,-S(2>.  The  rolu- 
uuyy  uuyy 

tions  for  the  two  modes  of  spin  vibrational  frequencies  are  then  found  from 
these  equations  to  be 


Ti2co2  = ( 12JS)2  (1  -R2/2)  (1  -y2) 

ri  a 

+ (12JS)2  (R2/2)  (1+  ya) 

- 24JSK(1  -R^/4)  (1  + ya)  . (3.10) 

We  now  make  use  again  of  the  low  temperature  approximation  for  y^  (i.e., 
equation  (2.8));  thus,  we  find  for  the  choice  of  the  negative  sign  and  the 
positive  sign  in  (3.  10)  that  the  frequencies  for  the  two  modes  of  State  II 
may  be  written  as 


:K2«^  = (12JS)4  (1  - 1 y 

l 4(12JS-K) 


(g(3H)‘ 


l(*  ■•) 


(3.  11) 


and 


■j.  2 2 

h w = 48JSK 


f , 12JS 

1'  K 

{■ 


2JS+  1}  __(.g^_H)  L 

^ A I 1 n TP  T.T  \ “ 


4(  1 2JS  - K) 


+ ( 1 2JS)  < 1+ 


(gft  H)‘ 


K _ (3+  K . 

12JS  TZJS  4(12js  -K)2 


^4tt  2 2 
3 ° 


(3-  12) 


t 


TR  198 


-17- 


respectively . It  may  be  readily  seen  that  the  forms  of  these  equations  are 
special  cases  of  the  more  general  form  of  equation  (3.  2).  Hence,  the  con- 
tribution of  the  mode  represented  by  (3.  1 1)  to  the  specific  heat  may  be  ex- 
3 

pressed  as  a T relationship  (similar  to  (3.4)),  while  the  specific  heat  due 

to  the  mode  represented  by  (3.  12)  must  be  expressed  in  a more  complicated 

form  (similar  to  (3.3)  with  "5  =0)*  For  a numerical  example  of  the  varia- 

tion  of  the  specific  heat  with  applied  field  at  a constant  temperature,  we  have 

chosen  the  conditions , 12JS/kT  = 15andK/12JS=  0.00111.  Thus.'E  = l,for 

o 

which  we  already  have  computed  the  specific  heat  for  State  I (as  indicated 
in  Figure  6b);  furthermore,  the  transition  between  States  I and  II  will  occur 
when  the  condition,  = 2[K(12JS  - K)]  ^/kT  = 0.999  is  satisfied.  Cal- 
culating the  specific  heat  contributions  of  the  two  modes  of  State  II,  repre- 
sented by  equations  (3.  II)  and  (3.  12),  we  find  that  with  increasing  applied 
field,  the  former  rises  very  slowly  while  the  latter  decreases  relatively 
quickly.  The  net  effect  is  a decrease  of  the  total  specific  heat  of  the  system. 
The  calculated  results  are  shown  in  Figure  8,  and  it  is  quite  evident  that 
there  is  a distinct  difference  between  States  I and  II  in  the  manner  in  which 
the  spin  wave  specific  heat  varies  with  a magnetic  field  applied  along  the 
anisotropy  axis.  We  should  also  note  that  at  the  point  of  transition  between 
States  1 and  II,  there  is  an  abrupt  discontinuity  in  the  total  specific  heat  of 
the  antiferromagnetic  system. 

Let  us  turn  now  to  the  situation  in  which  a magnetic  field  is  applied 
perpendicular  to  the  anisotropy  axis  (and  along  one  of  the  other  cube  edges) 
of  the  simple  cubic  antiferromagnetic  system.  In  this  case,  as  shown  in 
Appendix  B,  there  is  only  one  possible  static  equilibrium  spin  state  within 
the  applied  field  range  of  interest  (as  compared  to  the  two  for  the  case  just 
previously  discussed).  The  corresponding  spin  configuration  of  the  two 
sublattice  system  is  represented  schematically  in  Figure  7.  It  is  fairly 
obvious,  incidentally,  that  this  configuration  also  represents  the  only  stable 
Bpin  state  of  a simple  cubic  antiferromagnetic  structure  having  cubic  ani- 
sotropy, when  a small  field  is  applied  along  one  of  the  three  mutually  - 
perpendicular  anisotropy  axes.  With  reference  to  this  figure,  the  angle, 

9,  satisfies  the  condition 
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sin  0 = gf3H/2(12JS  + K)  , (3.13) 

and,  thus,  as  the  applied  field  is  increased  from  zero,  0 increases  mono- 
tonically  to  tt/2.  However,  for  reasons  previously  mentioned,  only  small 
values  of  0 may  be  expected  to  be  physically  attainable,  and  our  subsequent 
calculations  cf  the  spin  wave  specific  heat  will  be  restricted,  therefore, to 
this  region  of  applied  fields. 

The  starting  point  of  the  spin  wave  problem  is  once  again  the  equa- 
tion of  motion  (2.  1).  The  rectangular  components  of  the  spin  vectors  of 
the  two  sublattices  may  now  be  expressed  as 

S(l'=  j (Ssin0+S(1))  + j S(1)  +y  (Scos0  + S(1))  (3.14a) 

— Jx  x y y Jz  z 

S(2)=  j (S  sin0  + S(2))  + j S(2)  + j (-Scos  0 + S(2))  . (3.14b) 

The  applied  field  is  taken  to  be  in  the  x direction.  Again  assuming  that 
the  dynamic  spin  components  have  the  spatial-wave  form  of  (2.3),  we 
find  that  the  first  order  equations  of  motion  of  these  components  are 


12JSya  cos  6 + ( 12JS+2K)  cos  0 (3.  15a) 

W2^  - 12JSyacos  GS^  - ( 12JS+2K)  cos  0 S^2)  (3.15b) 

- 1 2 JSv  cos  OS*2*  - ( 12JS+2K)  cos  6 S(1) 
y la  x x 

+ 1 2JSy  sin  0S(2)  - [( 1 2 JS  + 2K)  sin  0 - g^H]  S(1)  (3.15c) 

QL  Z Z 

W2)=  12JSy  cos  0 s(1)  4-  ( 12JS+2K)  cos  0 S(2) 

y 1 a x x 

+ 1 2JSy  sin0S(1)  - [(12JS+2K)  sin  0 - gpH]  S(2)  (3.  15d) 

z z 

= -12JSya  sin  0S^  + [12JS  sin  0 - g£H]  (3.  15e) 

hS(z2)= -12JSyasin  + [12JS  sin  0 - gpH]  S^2)  . (3.  15f) 
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Since  the  magnitudes  of  the  spin  vectors  are  time  independent,  it  may  be 
deduced  from  equations  (3.  14)  that 


:(D 


X 

cos  0 


-s(1) 

z _ M) 

sin  9 v 


and 


S<2> 

X 

cos  ( 


s(2) 

z 

sin  8 


= S(2> 
v 


Making  use  of  (3.  13)  to  remove  0 from  equations  (3.  15),  we  find  that  in 
terms  of  the  y and  the  newly  defined  v components  of  spin,  equations  (3.  15) 


reduce  to  the  four  independent  equations, 

12JSYaS^2)  + (12JS+ 2K)  (3.16a) 

W2)=  -(12JS+  2K)S^2)  (3.16b) 

kS(1)=  - 12JS\  (1  -Q2/2)  S(2) 
y *a  H v 

- [12JS  + 2K(1 -Q2/4)]  S(1)  (3.16c) 

ri  v 

Wy2)=  12JSya(l  -Q^/2)S^1) 

+ [12JS  + 2K(1  -Q2/4)]S(2)  (3.1 6d) 

n v 


where  = gpH/(  12JS+K).  The  dynamic  spin  components  are  assumed  to 

vary  harmonically  with  time;  thus,  it  becomes  possible  to  combine  (3.  16a) 

with  (3.  16b),  and  (3.  16c)  with  (3.  16d)  to  form  one  pair  of  linear  equations 

in  S(1>-S<2>  and  S^+S(2)  and  an  other  in  S(1)+S(2*  and  S(1)-S(2^.  The 
vv  yy  vv  yy 

two  modes  of  spin  variational  frequencies,  corresponding  to  the  solutions 
of  these  two  pairs  of  equations,  are  found  to  be 

■k2w2  = 4K(  12JS  + K)  (1  - Q^/4) 

+ (12JS)2  (l-Q^/2)  (1  -y2) 

+ 12JS(  12JS  + K)  (Q2  /?)  (1+  y ) (3.17) 

which,  upon  use  of  (2.8)  as  the  low  temperature  approximation  for  yo  > is 
converted  to 
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ii2w4'  = 4K(12JS  + K)  I 1 — LsJLiil 

I 4(12JS  + K) 


+ (12JS)2(l  - 


4( 


jtfiai-l  (±£ ,/) 

1 2JS  + Ki  J \ / 


(3.18) 


and 


li2to  2 = 4K{  12JS  +K)J  - 1)  + 1 

L 4(  12JS  + K)^ 


(i2JS)2f  1 - (3  + K- i _umL.]C: 
L J 4(  1 2JS  + K)2  J \ 


4tt2  2N 
a 


, (3.19) 


corresponding  respectively  to  - and  + in  (3.  17).  In  the  absence  of  an  ap- 
plied field,  both  (3.  18)  and  (3.  19)  reduce  properly  to  (3.  2)  in  which  H 
has  been  set  to  zero.  The  two  modes  thus  become  degenerate.  From 
the  forms  of  (3.  18)  and  (3.  19),  it  follows  that  the  spin  wave  contribution 
of  each  mode  to  the  specific  heat  may  be  represented  by  an  expression 
that  is  similar  to  (3.  3)  with  ^ = 0.  For  comparison  with  our  previous 

calculations,  we  have  chosen  the  conditions,  12JS/kT  = 15  and  K/12JS  = 
0.00111,  for  which  to  work  out  a numerical  example  of  the  dependence 
of  the  specific  heat  (at  a fixed  temperature)  on  the  magnetic  field  applied 
perpendicular  to  the  anisotropy  axis.  The  calculated  results  are  shown 
in  Figure  9,  from  which  it  is  evident  that  for  an  increase  of  the  applied 
field,  the  rise  of  the  specific  heat  contribution  of  the  mode  represented 
by  (3.  18)  is  more  than  compensated  by  the  decrease  of  that  of  the  mode 
represented  by  (3.  19).  The  net  result  is  a decrease  of  the  total  specific 
heat. 


It  would  be  interesting  to  examine  together  the  effects  on  the  total 
spin  wave  specific  heat,  of  a magnetic  field  applied  parallel  to  the  ani- 
sotropy axis  and  of  a field  applied  perpendicular  to  this  axis.  Conse- 
quently, in  Figure  10,  we  have  superimposed  the  pertinent  curves,  which 
appear  separately  in  Figures  8 and  9,  for  the  conditions,  12JS/kT  - 15 
and  K/12JS  = 0.00111.  Moreover,  we  have  calculated  the  same  types  of 
curves  for  the  same  value  of  12JS/kT  but  for  K/12JS  = 0.00444.  These 


g/3H  /kT  - (WITH  T FIXED) 


g/3H  /KT  (WITH  T FIXED) 
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curves,  also,  are  shown  in  Figure  10,  It  becomes  evident  that  the  abrupt 
change  of  specific  heat  that  accompanies  the  transition  between  States  I 
and  II  (for  H parallel  to  the  K-axis)  is  greater,  the  larger  the  magnetic  ani- 
sotropy. Furthermore,  beyond  these  transitions  of  spin  state  (i.e.,  for 
higher  applied  fields),  a larger  anisotropy  causes  a bigger  difference  be- 
tween the  cpecific  heat  when  the  applied  field  is  parallel  to  the  anisotropy 
axis,  and  the  specific  heat  when  the  field  is  perpendicular  to  this  axis. 

If  the  magnetic  anisotropy  is  considered  to  be  negligibly  small.-  both 
equations  (3.  11)  and  (3.  18)  reduce  to 

lZw2=  [(12JS)2  - (gpH/2)2]  (4ir2a2/3)  , (3.20) 

and  both  (3.  12)  and  (3.  19)  become  equal  to 

k2u>2  = (gpH)2  + [( 1 2JS)2  - 3( gf3H/2)2]  (4n2a2/3)  . (3.21) 

Hence,  these  equations  represent  the  two  modes  of  spin  vibration  for  a 
magnetically  isotropic  simple  cubic  antiferromagnetic  system  with  a 
magnetic  field  applied  along  a cube  edge.  It  follows  from  previous  dis- 
cussion that  the  specific  heat  due  to  the  mode  represented  by  (3.20)  will 
3 

follow  a simple  T law,  while  the  temperature  dependence  of  the  specific 
heat  due  to  the  mode  represented  by  (3.21)  will  have  a more  complicated 

behavior  (described  by  an  equation  similar  to  (3,3)  with  B,j  = 0)*  Although 

n 

this  isotropic  case  may  not  have  much  physical  significance,  it  does  pro- 
vide an  interesting  extreme  case  for  which  the  field  dependence  of  the 
specific  heat  may  be  compared  to  that  of  the  system  with  a finite  magnetic 
anisotropy.  Therefore,  once  again  for  12JS/kT  = 15,  we  have  calculated 
the  specific  heat  contributions  of  the  modes  represented  by  ( 3.  20)  and  ( 3.  2 1), 
and  we  have  found  that  for  an  increase  of  the  applied  field  (at  a fixed  tem- 
perature), the  decrease  of  the  latter  more  than  compensates  for  the  in- 
crease of  the  former.  The  resultant  decrease  of  the  total  specific  heat 
is  quite  clear  from  Figure  10,  in  which  the  total  specific  heat  of  this  iso- 
tropic case  has  been  plotted  over  the  applied  field  range  of  interest. 
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4.  Comparison  with  Experiment 

Since  the  preceding  spin  w?.ve  analyses  were  developed  semiclarssi- 
cally  and  the  determinations  of  the  ground  state  spin  configurations  were 
purely  classical,  we  certainly  cannot  rely  rigorously  on  our  theoretical 
results  for  quantitative  comparison  with  experiment.  Nevertheles s , these 
results  could  be  taken  seriously,  were  it  possible  to  examine  them  with 
reference  to  any  thermal  or  magnetic  measurements  conducted  at  low 
temperatures  (where,  presumably,  the  spin  wave  theory  is  valid)  on  a 
ferrimagnetic  or  antiferromagnetic  material  having  a simple  cubic  struc- 
‘ ure . Unfortunately,  while  there  have  been  numerous  such  measurements 
on  ferromagnetic  metals  and  alloys,  only  a few  studies  have  been  made  of 
the  low  temperature  spin  wave  properties  of  a ferrimagnetic  or  antiferro- 
magnetic  structure  of  any  type  of  crystalline  symmetry.  This  is  particu- 
larly true  if  we  disregard  recent  spin  resonance  experiments,  which  are 
not  of  direct  interest  in  this  report.  Hence,  we  are  left  with  only  the 
possibility  of  making  some  qualitative  predictions  about  any  experiments 
that  are  suggested  by  our  simple  theory,  using  any  existing  and  pertinent 
experimental  information  for  starting  points  of  discussion. 

To  date,  the  main  experimental  evidence  for  the  existence  of  a 
3/2 

T law  has  been  the  temperature  dependence  of  the  saturation  magnetiza- 

17 

tier,  of  ferromagnetic  metals  at  low  temperatures.  It  is  significant, 
perhaps,  that  the  decrease  of  the  magnetization  of  gadolinium  was  found 

3/2 

to  be  proportional  to  T all  the  way  up  to  its  Curie  point,  while  the 

magnetizations  of  nickel  and  iron  at  intermediate  temperatures  follow 

2 1 8 
a T relationship  (which  is  explained  by  Stonei  on  a collective  electron 

theory  basis).  We  shall  refer  to  this  later. 

Of  several  papers  on  measurements  of  the  temperature  dependence 

19 

of  the  magnetization  of  various  ferrites,  the  most  recent  one  byPauthenet, 
containing  considerable  tabulated  data,  is  the  most  suitable  for  the  present 
discussion.  Unfortunately . except  for  a few  low  temperature  points  used 
to  establish  the  magnetizations  at  absolute  zero,  all  his  measurements  (on 
several  simple  ferrites  and  a family  of  mixed  ferrites)  were  made  above 
room  temperature.  Nevertheless,  we  shall  attempt  here  a crude  comparison 
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between  these  measurements  and  the  predictions  ol  our  simple  spin  wave 

theory.  Equation  (2.  18b)  is  the  relevant  expression,  if  the  effects  of 

magnetic  anisotropy  and  large  applied  fields  are  to  be  neglected.  The 

octohedral  (B)  and  the  tetrahedral  (A)  sites  in  the  actual  ferrite  struc- 
20 

tare  will  be  considered  roughly  equivalent  to  the  atom  sites  in  sub- 
lattices (1)  and  (2),  respectively,  of  our  simple  cubic  system;  thus,  ef- 
fectively, the  spins  of  the  atoms  of  the  two  B sites  of  each  unit  cell  will 
be  coupled  rigidly  together.  Hence,  for  the  simple  ferrites,  we  setS^ 

2.5  (i  . e.  , the  spin  of  one  ferric  ion),  and  = S^,  + 2.5,  where  S^,  is  the 
experimentally  determined  average  magnetic  moment  per  atom  divided 
by  2p  . For  the  mixed  nickel -zinc  ferrites  (i  . e .,  (NiO)  ^ ^(ZnOJ^Fe^O , 
we  shall  allow  for  the  apparent  substitution  of  non-magnetic  zinc  ions 
for  ferric  ions  in  the  A sites  and  the  transfer  of  these  ferric  ions  (with 
reversed  spin)  to  the  B sites , by  setting  = 2.5  ( 1 -x)  and  S^=S,p  + 2.5(l-x). 
This  can  be  only  even  approximately  valid  for  x<  0.5  since  the  magnetic 
moment  per  ferrite  molecule  actually  decreases  for  larger  values  of  x. 

Let  us  now  find  a relationship  between  the  exchange  integral  in 

(2.  18b)  and  the  Curie  temperature,  T , According  to  the  modern  version 

21  C 

of  the  Weiss  field  theory,  we  may  write  for  the  z components  (i.e., 
those  in  the  direction  of  the  applied  field)  of  the  magnetizations  of  the 
two  sublattices: 

M(1)=  NgpS(1)=  NgPS,  • Bc  (12JSrS(Z)/kT) 

Z Z 1 £>  J i Z 

M(Z)=  Ng(3S(Z)=  NgPS,  • Bc  (12JS_S(1)/kT)  . (4.1) 

z z l c>  z 

Remembering  that  the  Brillouin  function,  Bg(u)  — (S+l)u/3S  for  u <<C  1 , 

we  obtain  from  (4. 1)  two  linear  equations  in  and  S^2^,  whose  com- 

z z 

patibility  condition  fixes  Tc  as  follows: 

kTc  = 4J  [SjS^Sj+lMS^l)]  1/2  . 

Consequently,  (2.  18b)  may  be  rewritten  * s 


(4.2) 
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A .M/M  = A.  .(TAT  )3/2 
o o M c 

Am  = 0.  ll7(S1-S2)1/2L(S1+i)(S2+l)/S1S2]3/4  (4.3) 

when  anisotropy  and  applied  field  effects  are  neglected. 

Starting  with  Pauthenet's  values  for  Tc  and  S^,,  we  have  calculated 
S.,  S2>  J,  and  for  cur  simple  cubic  models  for  the  simple  ferrites  of 
nickel,  cobalt,  iron,  and  manganese,  and  for  various  mixed  nickel-zinc 
ferrites.  The  results  are  listed  in  Table  1.  Furthermore,  we  have  con- 
verted Pauthenet's  data  into  equivalent  Aj^  versus  T/Tc  points,  whichare 
shown  in  Figure  11.  The  solid  curves  through  these  points  have  been 
restricted  to  the  temperature  range  over  which  we  believe  the  experimental 
points  are  dependable.  The  values  of  Aj^  which  we  previously  calculated 
and  listed  in  Table  1,  are  represented  in  Figure  11  by  the  horizontal  long- 
dashed  lines . 

By  a different  manipulation  of  his  data,  Pauthenet  was  able  to  con- 
clude that  A M/M  at  low  temperatures  for  nickel,  cobalt,  and  iron  ferrites 

°Z  ° 3/2 

followed  a T law,  while  that  of  manganese  ferrite  varied  as  T , The 

short  dashed  lines  in  Figure  11  correspond  quantitatively  to  these  tempera- 
ture variations.  However,  from  this  figure,  it  appears  to  us  that  for  all 
four  simple  ferrites,  as  well  as  for  the  mixed  ferrites,  it  is  quite  likely 
that  the  limiting  values  of  A..,  as  T — >0,  are  other  than  zero.  If  they 

iVl 

should  approach  their  limiting  values  along  a horizontal  tangent  in  Figure 
3/2 

11,  the  T law  for  low  temperatures  would  be  confirmed.  Moreover,  as 

in  the  cases  of  metallic  iron  and  nickel,  the  temperature  dependence  of 

AqM/Mq  for  nickel,  cobalt  and  iron  ferrites  probably  shifts  gradually  from 

T^/2  to  T^,  as  the  temperature  is  increased  from  a very  low  value.  The 

relative  magnetization  of  manganese  ferrite,  however,  like  that  of  metallic 

3/2 

gadolinium,  probably  follows  a T law  over  the  entire  ferrimagnetic 
temperature  range.  Our  calculated  low  temperature  values  of  Aj^  are  all 
undoubtedly  too  low.  However,  particularly  for  the  mixed  ferrites,  their 
relative  magnitudes  appear  to  be  in  the  right  sequence.  This  is  as  much 
as  could  be  expected  in  view  of  the  drastic  approximation  made  in  fitting 
the  actual  ferrites  into  "equivalent"  simple  cubic  models.  A spin  wave 


L 
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Table  I 


YO  • 


S?  Tc(°K)  J/k 


Ni  1.11  3.61  2.50  870  18.0  0.190 

Co  1.97  4.47  2.50  769  13.1  0.245 

Fe  2.04  4.54  2.50  847  14.3  0.250 


- 3 -1 
q T 1 millijoules  2 


Ml 


Mr  2.20  4.70  2.50  576 


9.4  0.258 


K - mole 


0 132 
0.366 
0.332 
0.660 


(Ni°)  i _x(  Zn0)x'  Fe203 


X S, 


S2  Tc(uK)  J/k 


c.  T ^ millijoules  ^ Oj 


K - mole 


0.1 

1.51 

3.76 

2.25 

814 

•»  •n  n 
11.0 

0.226 

0.250 

0.2 

1 ,89 

3.89 

2.00 

755 

17.7 

0.259 

0.378 

0.3 

2.24 

3.99 

1.75 

671 

17. 1 

0.290 

0.600 

0.4 

2.48 

3.98 

1.50 

618 

17.9 

0.320 

0.830 

0 5 

2.72 

3.97 

1.25 

548 

18.4 

0.355 

1.206 
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analysis  more  appropriate  for  the  spinel  structures  of  the  ferrites , should 
make  more  dependable  calculations  possible;  we  propose  to  study  this 
problem  in  the  near  future. 

Needless  to  say,  accurate  low  temperature  measurements  of  the 
saturation  magnetizations  of  the  various  ferrites  are  badly  needed.  The 
accuracy  required  to  determine  the  temperature  dependence  of  un- 

ambiguously, however,  is  extremely  difficult  to  attain  at  low  temperatures 
where  a very  small  change  of  a very  large  quantity  must  be  measured. 

Here,  we  suspect,  lies  the  main  advantage  of  low  temperature  measure- 
ments of  the  spin  wave  contribution  to  the  specific  heat.  Another  advantage 
of  specific  heat  measurements  over  those  of  magnetization  is  the  fact  that 
any  appreciable  magnetic  anisotropy  of  the  ferrimagnetic  structure  maybe 
expected  to  have  a greater  effect  on  the  temperature  variation  of  magnetiza- 
tion than  on  that  of  the  specific  heat  (sea  Figure  3).  The  large  applied 
fields  required  for  the  magnetic  measurements  would  tend  to  aggravate 
this  situation.  Hence,  the  experimental  results  for  the  temperature  de- 
pendence of  the  specific  heat  can  probably  be  more  easily  interpreted  in 

3/2 

terms  of  the  coefficient  of  a simple  T relationship  (see  equation  (2.  15b)), 
than  can  the  results  for  the  temperature  dependence  of  magnetization. 

In  many  ferromagnetic  metals  and  alloys  whose  low  temperature 

specific  heats  have  been  measured,  the  electronic  specific  heat  due  to 

the  energy  band  structure  of  the  cl  electrons,  has  been  large  enough  to 

mask  any  possible  spin  wave  contribution,  We  have  found,  however,  that 

the  spin  wave  specific  heat  of  a ferrimagnetic  simple  cubic  structure  also 
3/2 

follows  a T law  (see  equations  (2.  15))  and  is  thus  experimentally  dis- 
tinguishable from  that  due  to  lattice  vibrations  (which,  of  course,  ic  pro- 
3 

portional  to  T at  low  temperatures).  It  seems  reasonable  to  assume  that 
this  is  true  for  ferrimagnetic.  materials  of  any  crystalline  structure. 

Hence,  the  ferrites,  whose  electronic  specific  heats  are  probably  very 
small,  appear  to  be  the  obvious  materials  for  spin  wave  specific  heat 
measurements . 

Using  (2.  15b)  and  (4.  2),  we  have  calculated  the  spin  wave  specific 
heats  for  the  different  ferrites  discussed  previously.  The  resulting 
specific  heat  values,  which  have  been  listed  in  Table  1,  are  per  mole, 
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where  the  molecule  has  been  taken  co  consist  of  two  adjacent  lattice  sites 

of  the  simple  cubic  model.  Note  the  comparatively  high  spin  wave  specific 

heats  of  the  mixed  ferrites  that  have  relatively  low  Curie  temperatures 

(but  whose  exchange  integrals,  incidentally,  are  not  very  different  from 

those  of  the  simple  ferrites).  However,  it  is  the  magnitude  of  the  spin 

wave  specific  heat  relative  to  that  of  the  lattice  specific,  heat,  that  is  of 

prime  experimental  importance.  About  the  only  specific  heat  measure- 

22 

ments  on  any  of  the  ferrites,  to  date,  have  been  those  of  Millar  on  iron 
ferrite  (magnetite)  down  to  liquid  air  temperatures.  Assuming  that  at 
very  low  temperatures,  the  Debye  temperature  of  iron  ferrite  is  about 
the  same  (i.e.,  about  565°K)  as  at  higher  temperatures,  we  find  that  the 
spin  wave  and  lattice  vibration  contributions  to  the  specific  heat  are  both 
equal  at  2.7  K to  about  1.5  millijoules/  K mol.  At  temperatures  higher 
than  2.7°K,  the  lattice  specific  heat  will  predominate.  Hence,  the  specific 
heat  measurements  should  be  carried  out  down  to  liquid  helium  tempera- 
tures. The  low  specific  heat  will  necessitate  extreme  care  in  the 
calorimetry.  An  interesting  supplementary  experiment  at  these  tempera- 
tures, preferably  on  a single  crystal  ferrite  specimen,  would  be  the 
measurement  of  the  dependence  of  the  spin  wave  specific  heat  on  the 
strength  of  an  applied  magnetic  field. 

If  the  effects  of  anisotropy  and  applied  fields  are  ignored,  the  spin 

wave  specific  heat  of  a simple  cubic  antiferromagnetic  structure  has 

3 

been  shown  in  Section  3 to  be  proportional  to  T . Thus,  it  would  not  be 
easily  separable  experimentally  from  the  lattice  specific  heat  at  low 
temperatures,  unless  the  spin  wave  contribution  were  very  large  due  to 
the  small  magnitude  of  the  exchange  interaction  integral  (see  equation 
(3.4)).  The  Curie  or  Neel  temperature  for  the  antifer  rooragnetic  struc- 
ture may  be  expressed  by  (4.  2)  if  we  set  = S;  the  equation  then 

becomes  simply 

kTc  = 4JS(S+ 1)  . (4.4) 

Hence,  it  follows  that  spin  wsvi  specific  heat  measurements  should  be 
possible  with  artifer romagnetic  materials  having  low  Curie  temperatures 
(and  low  electronic  specific  heats). 


r 
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Specific  heat  measurements  were  recently  completedby  Friedberg 

on  cupric  chloride.  Up  to  the  Curie  temperature  (about  4.  3°K),  the  total 

specific  heat  of  this  antiferromagnetic  compound  was  found  to  be  over  one 

hundred  times  as  large  as  the  estimated  lattice  specific  heat  and  was 

3 

roughly  proportional  to  T . However,  since  these  measurements  were 
made  down  to  only  1.5°K,  it  is  almost  certain  that  the  anomolously  large 
specific  heat,  as  well  as  its  cubic  temperature  dependence,  are  associated 
primarily  with  an  order -disorder  transformation  of  the  spin  state,  rather 
than  with  small  perturbations  of  an  ordered  spin  state.  Hence,  to  deter- 
mine the  normal  spin  wave  specific  heat  of  cup.i ic  chloride,  it  would  be 
necessary  to  extend  these  measurements  down  to  temperatures  in  the 
region  of  T^/10  (i.e.  , about  0.5°K)  and  study  the  variation  of  the  specific 
heat  over  a very  restricted  temperature  range. 

Although  most  other  antiferromagnetic  salts  have  higher  Curie  points, 
we  find,  when  we  substitute  (4.4)  in  (3.4)  and  obtain 

C = 0.507  k(S+ l)3  (T/T  )3  , (4.5) 

ii  C 

that  a larger  electronic  spin  of  the  magnetic  ions  would  tend  to  counter- 
balance the  effect  of  a higher  Curie  point  on  the  magnitude  of  the  spin 
wave  specific,  heat.  Let  us,  therefore,  make  a rough  estimate  of  the  spin 
wave  specific  heat  of  ferrous  chloride,  whose  Curie  point  is  about  24°K, 
by  assuming  that  the  spin  associated  with  each  lattice  site  of  our  simple 
cubic  model  is  2 (i . e . , the  spin  of  one  ferrous  ion).  We  substitute  this 

value  of  S into  equation  (4.5)  and  obtain  for  the  spin  wave  specific  heat: 

3 o 4 

8.2  T millijoules/(  K)  mole,  (a  molecule  is  considered  equivalent  to 

one  lattice  site  of  the  simple  cubic  model).  Let  us  then  suppose  that  the 

lattice  specific  heat  of  ferrous  chloride  is  not  very  different  from  that 

3 o 4 _ 

of  cupric  chloride  (i.e.,  about  1.1  T"'  millijoules/(  K)  mole).  Hence, 

the  spin  wave  specific  heat,  due  to  its  large  relative  magnitude,  should 

be  quite  separable  from  the  lattice  specific  heat.  Furthermore,  we  find 

from  (4.4)  that  J = k,  and  if  we  assume  that  K— 0.  1 k,  it  follows  that 
1/2 

2[K(12JS  + K)]  — 3.  1 k.  Thus,  with  reference  to  Figures  5a  and  5b, 

it  appears  that  for  the  temperatures  at  which  the  specific  heat  measure- 
ments on  ferrous  chloride  should  be  carried  out  (i.e.,  in  the  region  of 
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T / 10  or  2.4°K),  the  spin  wave  specific  heat  should  be  well  out  of  the  ex- 
ponential temperature  dependence  region. 

24 

Magnetic  measurements  on  single  crystals  of  manganese  flouride 
25 

and  cupric  chloride  have  confirmed  Neel's  prediction  that  a transition 
of  spin  state  occurs  when  a magnetic  field  applied  along  the  anisotropy  axis 
of  an  antiferromagnetic  crystal,  attains  a certain  critical  value.  It  ■would 
be  most  interesting,  therefore,  to  make  the  accompanying  specific  heat 
measurements  and  learn  if  the  magnetic  field  dependence  of  the  spin  wave 
contribution  can  be  represented  by  a curve  similar  to  that  in  Figure  8. 

For  the  sake  of  comparison,  these  specific  heat  measurements  should 
also  be  made  with  the  field  applied  perpendicular  to  the  anisotropy  axis 
(see  Figure  10). 
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Appendix  A 

The  Effects  of  Next-to-Neare at  Neighbor  Exchange  Interactions 


In  the  simple  cubic  atomic  spin  structure  shown  in  Figure  1,  each  atom 
has  twelve  next-to-nearest  neighbors,  all  of  which  lie  on  the  same  spin  sub- 

3 

lattice.  To  be  consistent  with  the  cur  rent  theory  on  the  properties  of  ferrites, 
we  assume  that  the  next-to-neare  st  neighbor  exchange  interactions  for  both 
sublattices  are  negative,  similar  to  but  considerably  weaker  than  the  nearest 
neighbor  interactions.  Hence,  the  net  electronic  spins  of  nearest  neighbor 
atoms  are  still  mutually  antiparallel,  although  the  net  exchange  force  on  each 
spin  is  somewhat  smaller  than  it  would  be  in  the  absence  of  next-to-nearest 
neighbor  interactions.  We  may  expect,  therefore,  that  the  introduction  of  nega- 
tive next-to-nearest  neighbor  interactions  into  our  analysis  will  result  in  an 
increase  of  the  spin  wave  specific  heat  of  the  system.  It  may  be  of  interest, 
however,  to  examine  their  effects  more  quantitatively. 


We  let  -J^  and  -J^  be  the  exchange  integrals  for  the  next-to-nearest  neighbor 
interactions  between  spins  on  sublattices  (1)  and  (2),  respectively.  As  before, 

-J  is  the  exchange  integralfor  the  nearest  neighbor  interactions . Thus,  J^,  J^, 
and  J are  all  positive  quantities.  Furthermore,  J>>Jj  or  J^.  If  we  now  ->.ssume 
that  the  equation  of  motion  of  a particular  spin  vector  has  the  form  of  (2.1), 
and  then  use  (2.2)  and  (2.3),  we  find  that  the  equations  of  motion  for  the  spin  com- 
ponents are  the  same  as  (2.4)  except  that  to  (a),  (b),  (c),  and  (d)  must  be  added: 

-241,3.(1-7  ' )S(1),  +24J  S,(l-  y 1 )S(2),  +24J.S,(1-  y ' )S(1),  and  - 24J„S,(1-  y ' )S(2), 

1 r a y 2'  a y 11'  a x 2 2 ax’ 


respectively,  where  y^  = (cos  2Tfaycos  ^ira^  +cos  2irayCOS  2iraz+cos2irazcos 
When  (2.5a)  is  used,  this  is  equivalent  to  adding  -24J^S^(1-  y^)S^  and 


+24J2S2<1^d)SV 

gx  : 8 2 = g and  K! 


to  (2.  6a)  and  (2.  6b),  respectively.  Thus,  again  setting 
= K2  = K-  we  obtain  for  the  spin  wave  frequencies: 


-29- 
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tw  = gpH  - 6J(S.  - S2)  - 12(J1S1  - J2S2)(1  - 

+ - 5?)  + 12{J1S1  - J2S2)(1  -y^)}  2 + {12J)2S1S2{1  -yf) 

+ 2K  {l2J(S1  + S2)  + 2K  - 24(J1S1  + - Y^)] 

- 24{J1S^  + J2S^)(12J)(1  -y'J 


2l1/2 

+ (24J1S1){24J2S2){1  - Y;)  I - 


(A.  1) 


For  spin  waves  of  very  long  wavelength,  a is  very  small,  and,  therefore, 

2 2 2 
1 - Y1  — 4tt  q /3,  which  is  also  the  small -a  value  of  1 - -y 
■ a ’q 

In  the  ferrimagnetic  case,  6J(S^  - S2)  may  be  expected  to  be  by  far  the 
the  dominant  term  in  {A.  1);  hence,  the  plus  sign  must  be  chosen  for  the 
"acoustical"  mode.  Furthermore,  we  can  now  expand  the  square  root  in 
(A.  1)  binomially,  and,  neglecting  higher  order  terms  in  a,  we  get 


fS  + S?  + K/bJ  ' 

U«gPH  + 2K<  } 

L 1 2 J (A.  2) 

+ [JS!S2  - 2J1S1  K +£}  - 2J2S2  (S2  + &]]  hA2/3>. 


A comparison  of  (A.  2)  with  iZ.  9)  reveals  that  the  only  effect  of  next-to- 

nearest  neighbors  interactions  on  the  expression  for  w is  a decrease  of  the 

coefficient  of  a2.  It  can  be  easily  shown  that,  correspondingly,  both  the 

specific  heat  and  the  temperature  variation  of  the  magnetization  are  increased 

2 

by  a factor  equal  to  the  3/2  power  of  the  ratio  of  the  coefficient  of  a in  (2.9) 
to  that  in  (A. 2).  The  dependences  of  Ca  and  AoM/Mq  on  the  applied  magnetic 
field,  however,  are  still  described  by  (2.15a)  and  (2.18a),  respectively. 

For  the  antiferrornagnetic  case,  we  set  = S and  = J2;  as  a 

result,  we  find  from  (A.l)  that  the  expression  analogous  to  ('2.2)  may  be 
written  as 


t 


r 
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^i|c 


4J, 


4K(12JS  +K)  +( 12JS)  1 - 


/,  . K\ 

^ +5J Sj 


/4tt^ 

V3 


-A 


1/2 


J-rrflU 


(A.3) 


Jt  follows  that  the  introduction  of  next-to-neare st  neighbor  exchange  inter- 
actions causes  an.  increase  of  the  specific  heat  of  the  antiferromagnetic 
system  by  the  factor  - 4J^{1  + K/6JS)/jj  As  in  the  ferrimagnetic 

case,  the  field  dependence  of  the  specif1''  is  not  changed  by  these  addi- 

tional interactions . 


Since  all  next-to-neare  3t  neighbor  atoms  in  our  simple  cubic  model  lie 
in  the  same  sublattice,  their  electronic  spin  vectors  have  parallel  ground- 
state  orientations  whether  the  system  is  in  State  I or  in  State  II  (see  Figure  7). 
Hence,  the  criterion  for  the  transition  of  the  antiferromagnetic  system  from 
one  spin  state  to  the  other  is  net  affected  by  next-to-nearest  neighbor  inter- 
actions . 


L 
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Appendix  B 

The  Static  Equilibrium  States  of  the  Antiferromagnetic  Spin  System 

It  will  be  assumed  that  the  electronic  spin  vectors  of  the  simple  cubic 
atomic  system  are  subjected  only  to  forces  arising  from  negative  exchange 
interactions  between  nearest  neighbor  spins  and  from  interactions  with  a 
magnetic  amsotropy  field  and  a uniform  applied  magnetic  field.  The  corre- 
sponding energies  will  be  expressed  as  the  purely  classical  equivalents  of 
the  energy  terms  in  equation  (2.1).  The  g-factors  as  well  as  the  anisotropy 
constants  associated  with  the  two  sublattices  will  be  considered  equal. 

We  are  interested  in  determining  the  static  equilibrium  positions  of  the 
spin  vectors,  first,  when  the  applied  field  is  parallel  to  the  anisotropy  axis. 

It  is  fairly  obvious  in  this  case  that  the  spin  vectors  must  lie  in  planes  that 
are  parallel  to  each  other  and  to  the  anisotropy  axis.  The  problem,  there- 
fore, may  be  formulated  with  reference  to  the  schematic  two-dimensional 
vector  diagram  of  Figure  12,  in  which  the  orientations  of  two  neare st  neighbor 
spin  vectors  are  defined  by  thi£i  ang  les,  "X  and  0.  The  energy  of  the  system 
of  N atomi  may  be  written  then  as 

E = - 6NJS2cos  29  - -^NgpHS[  cos(X  - w)  + cos(ir  - X - 0)] 

+ -^NKS[  sin2(X  - 0)  + sin2(Tr  - X - 0)] 
or  E/NS  = 6JS(2sin20-l)  - g(3H  sin  0 sin  X 

+ K(sin20  + sin2  X - 2 sin2©  sin2X).  (B.l) 

It  is  easily  established  that  as  long  as  K<  6JS  (which  is  certain  to  be 
satisfied  in  practice),  the  value  of  0 that  minimizes  (B.l)  may  be  expressed  by 

8ine  = (g£H/2)  sin  X > (B.2) 

12JS  + K(  1 - 2 sin  X) 

for  which  (B  1)  becomes 

E/NS  = -6JS  + K sin2X (gpH/2)  sin  X (B.3) 

12JS  + K(1  ■■  2 sin  X) 


FIGURE  12 
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It  is  then  found  that  with  respect  to  "X  , (B.  3)  has  stationary  values  for 

sin  X = 0,  sin  9 = 0;  (B.4a) 

sin  X = 1,  sinG  = gpH/2(12JS  - K);  (B.4b) 


and 


2 

s in  X 


1 2 JS  + K 
2K 


1 - 


(gpH/2) 


[ K(  1 2JS  + K)] 


Vz 


1 2 JS  + K . 2q 
— sin  9, 


K 


(B . 4c) 


for  which  the  corresponding  energies  are  respectively, 

E/NS  = - 6JS,  (B.  5a) 

E/NS  = - 6 JS  + K - (g(3H/2)Z/<12JS  - K),  (B.  5b) 


and  E/NS  = - 6 JS  + — tK 


(flgH/2), 

[ K(  1 2JS  + K)]  1/2 


(B.  5c) 


Equation  (B.  5c)  clearly  represents  a maximum  energy  (which,  in  view  of 
(B.4c),  can  occur  only  within  a restricted  range  of  H).  Hence,  to  obtain 
the  absolute  minimum  energy,  we  must  compare  (B.  5a)  and  (B.  5b).  It 
follows  immediately  that 

1/?. 

gpH  = 2[K(12JS  - K)]  (B.  6) 

defines  a critical  value  of  H,  below  which  State  1 (defined  by  (B.  4a)  and 
(B.  5a))  will  be  the  lower  energy  state,  and  above  which  State  II  (defined 
hy  (B.  4b)  and  (B.  5b))  will  have  the  lower  energy. 

In  Figure  12,  the  energy  (as  expressed  in  (B.3))  has  been  plotted  as 
a function  of  X.  for  different  values  of  H;  the  unrealistically  high  value  of 
0.25  was  chosen  for  K/12JS  in  order  that  the  characteristics  of  equation 
(B.3;  that  are  dependent  on  K,  be  brought  out  more  distinctly.  It  is  clear 
from  this  figure  that  for  values  of  applied  field  other  than  those  close  to 
the  condition  for  spin  state  transition,  both  States  I and  II  may  be  expected 
to  be  reasonably  stable.  However,  when  the  field  is  very  near  the  critical 
value  defined  by  (B.6),  States  I and  II  are  favored  about  equaiiy  and  are 
separated  by  a low  energy  barrier  (see  Figure  12).  An  expression  for  the 
height  of  this  barrier  is  easily  determined  by  the  substitution  of  (B.6)  into 
(B.5c);  it  is  thus  found  that 


r 


I 
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( AE) 

— rn; — 

1ND 


barrier  = 


12JS  - [(12JS)2 


o'-  K2/24JS  forK<^12JS. 


/ »»  n\ 
[13.  I / 


In  most  physical  cases  of  antiferromagnetic  crystals,  this  energy  barrier 
is  so  low  that  a small  thermal  energy  is  sufficient  to  overcome  it.  Conse- 
quently, transition  between  States  I and  II  should  take  place  very  readily 
when  the  applied  field  is  equal  to  the  critical  value  defined  by  (B.6). 

Let  us  now  find  the  static  equilibrium  positions  of  the  spin  vectors  when 
the  field  is  applied  perpendicular  to  the  anisotropy  axis.  In  this  case,  al- 
though the  spin  vectors  must  again  lie  in  parallel  planes  due  to  the  exchange 
forces  between  them,  it  is  not  immediately  obvious  that  these  planes  must 
be  parallel  to  both  the  anisotropy  axis  and  the  axis  of  the  applied  field.  Hence, 
in  a plane  perpendicular  to  H,  let  us  turn  the  K-axis  through  an  arbitrary  angle, 
6,  from  its  position  shown  in  the  spin  vector  diagram  of  Figure  13.  With 
reference  to  this  schematic  diagram,  the  positions  of  the  spin  vectors  in  the 
plane  of  the  applied  field,  are  defined  by  the  angles,  X and  0.  The  energy  of 
this  system  rnay  then  be  written  as 


E = - 6NJS2  cos  29  - iNgpHS[  cofi(  X - 0)  + cos(tt-X  -0) 

+ -^NKStsin2^  ^ + sin2^^) 


where  cos  4^  = cos^cos(^ -X  + G) , and  cosi^  - cos<4cos(^  -X-  0). 
or  E/NS  = 6JS(2sin20  - 1)  - g(3Hsin9sinX 


K[1  - 


2,.  ■ 2n 

cos  p(sin  0 


. 2-* 
sin  a.  - 


2 sin^G  sin2X) • 


(B.8) 


It  follows  from  this  last  expression,  that  the  energy  is  minimized  with  respect 
to  when  6 is  set  to  zero.  Hence,  the  two-dimensional  schematic  diagram  of 
Figure  13  describes  the  system  adequately  for  the  remaining  discussion.  We 
proceed  to  find  that  (B.8)  is  minimized  with  lespect  to  0 when 


(gftH/  2) sin  X 


12JS  - K(1  - 2 sin  X) 


sin  0 = 


(B.9) 
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ln  whichcase,  (B.8)  becomes 

2 2 

E/NS  = - 6JS  + Kcos2 3 4X (B.10) 

1 2JS  - K(  1 - 2 sin  X) 

Finally,  by  differentiating  (B.10)  with  respect  toX,  we  find  that  as  long  as 
K<  12JS,  the  energy  is  minimized  by  X = it/ 2,  for  any  value  of  applied  field. 
Thus,  (B.9)  and  (B.10)  are  reduced,  respectively,  to 

sin©  = gj3H/ 2(12  JS  + K)  (B.ll) 

and  E/NS  = - 6JS  - (g|3H/2)2/(  1 2JS  + K).  (B.12) 

For  comparison  with  Figure  12,  the  energy  (as  expressed  in  (B.10)) 
has  been  plotted  in  Figure  13  as  a function  of  X for  K/  12JS  = 0.25  and 
various  values  of  H.  Clearly,  when  the  applied  field  is  perpendicular  to 
the  anisotropy  axis,  only  one  spin  state  of  the  system  is  possible,  and  it 
is  a stable  static  equilibrium  state  for  any  value  of  the  applied  magnetic 
field. 

Our  determination  of  the  spin  equilibrium  states  for  the  two  directions 
of  the  applied  field  relative  to  the  anisotropy  axis,  is  similar  but  more 
exact  than  Neel's  original  analysis.  ^ Recently,  Gorter  and  Haantjes2^* 
in  their  study  of  this  problem,  postulated  an  anisotropic  exchange  inter- 
action between  adjacent  spins  and  thus  included  the  magnetic  anisotropy  energy 
implicitly  in  their  expression  for  the  exchange  energy.  Their  results, 
nevertheless,  are  qualitatively  similar  to  ours. 

References 

1.  F.  Keffer  and  C,  Kittel,  Phys . Rev.  8_5,  329  (1952);  J.  Ubbink,  Physica 
19,  9 (1953);  -R.  K.  Wangsness,  Phys.  Rev.  9_3,  68(1954);  and 
references  therein. 

2.  P.  W Anderson,  Phys.  Rev.  694(1952);  R.  Kubo,  Phys.  Rev. 

87 . 568(1952);  J.  M.  Ziman,  Proc.  Phys.  Soc.  A 65,  549  and  548 
(19  52);  and  references  therein. 

3.  L.  Neel,  Ann.  Physique  3^  137  (1948). 

4.  J.  H.  Van Vleck,  Phys.  Rev  52,  1178(1937). 


TR198 


-36- 


References 


5 C.  Hearing  and  C Kittel,  Phys.  Rev.  (H,  869(1951). 

6.  R.  K Wangsness,  Phys.  P.ev.  9^ , 1085(1953). 

7.  F.  Seitz,  Modern  Theory  of  Solids  (McGraw  Hill,  New  York.,  1940), 
chap.  3,  sec.  18. 

8.  H.  B.  Dwigh*,  Tables  of  Integrals  and  Other  Mathematical  Data 
(Macmillan,  New  York,  1947),  p.  201 . 

9.  C.  Miller,  Z.  Physik  82,  559(1933). 

10  F.  Bloch,  Z.  Physik  6_1_,  206  (1930). 

11.  N.  F.  Mott  and  H.  Jones,  Theory  of  the  Properties  of  Metals  and 
Alloys  (Oxford  University  Press,  1936) , pp.  236-239. 

12.  T.  Holstein  and  H.  Primakoff,  Phys.  Rev.  58,  1098  (1940). 

13.  C.  Kittel,  Phys.  Rev.  82,565(1951). 

14.  Keffer,  Kaplan,  and  Yafet,  Am.  J.  Phys.  _21_,  250(1953). 

15.  H.  B.  Dwight,  op.  cit.  , p.  200. 

16.  L.  Neel,  Ann.  Physique^,  232(1936). 

17.  M.  Fallot,  Ann.  Physique  305(1936);  Elliott,  Legvold,  and  Spedding, 
Phys.  Rev.  9J_,  28(19  53). 

18.  E.  C.  Stoner,  Proc.  Roy.  Soc . A 165,  372  (1938). 

19.  C.  Guillaud  and  M.  Roux,  C.  R.  Acad.  Sci.  229,  1133  (1949); 

C.  Guillaud  and  H.  Creveaux,  C.  R.  Acad.  Sci.  230,  1256  and 
1458  (1950);  R.  Pauthenet,  Ann.  Physique  1_,  710_0'i52). 

20.  E.  J.  W.  Verwey  and  E.  L.  Heilmann,  J.  Chem.  Phys.  15,  174(1947); 
Verwey,  Haayman,  and  Romeijn,  J.  Chem.  Phy.ii.  1 5,  181X1947);  Shull, 
Wollan,  and  Koehler,  Phys.  Rev.  84,  912(1951). 

21.  C.  Kittel.  Introduction  to  Solid  State  Physics  (Wiley.  New  York.  1953), 
chaps  9 and  1 0 . 

22.  R.  Millar,  J.  Am.  Chem.  Soc.  50,  1875(1928)  and  51 , 215(1929). 

23.  S.  A.  Friedberg,  Physica  1(5,  714(1952). 

24  J.  W.  Stout  and  M.  Griffel,  J.  Chem.  Ph>  . 1455(1950). 

25.  van  den  Handel,  Gijsman,  and  Poulis,  Physica  18,862(1952). 

26.  C.  J.  Gorier  and  J.  Haantjes,  Physica  18,  285(1952). 


I 

Technical  P.eport a 

DISTRIBUTION  LIST 

2 

Chief  of  Naval  Research  (4  27) 
Department  of  the  Navy 
Washington  25,  D.  C. 

1 

Chief  of  Naval  Re search(460) 
Department  of  the  Navy 
Washington  25,  D.  C. 

1 

Chief  of  Naval  Research  ^421) 
Department  of  the  Navy 
Washington  25,  D.  C. 

6 

Director  (Code  2000) 

Naval  Research  Laboratory 
Washington  25,  D.  C. 

2 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 
150  Causeway  Street 
Boston,  Massachusetts 

1 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 

1000  Geary  Street 

San  Francisco  9,  California 

1 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 
1030  E.  Green  Street 
Pasadena,  California 

1 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 
The  John  Crerar  Library  Building 
86  East  Randolph  Street 
Chicago  1,  Illinois 

1 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 

346  Broadway 

New  York  13,  New  York 

3 

Officer  -in -Char  ge 
Office  of  Naval  Research 
Navy  No.  100 
Fleet  Post  Office 
New  Y ork,  N.  Y . 

-i- 

L 



) 


1 


l 


l 


2 


1 


1 


1 


1 


1 


2 


1 


1 


50 


Chief,  Bureau  of  Ordnance  (Re4) 

Navy  Department 
Washington  25.  D.  C. 

Chief,  Bureau  of  Ordnance  (AD-3) 

Navy  Department 
Washington  25,  D.  C. 

Chief,  Bureau  of  Aeronautics  (EL-i) 
Navy  Department 
Washington  25,  D.  C. 

Chief,  Bureau  of  Ships  (810) 

Navy  Department 
Washington  25,  D.  C. 

Chief  of  Naval  Operations  (Op-413) 
Navy  Department 
Washington  25,  D.  C. 

Chief  of  Naval  Operations  (Op- 20) 

Navy  Department 
Washington  25,  D.  C. 

Chief  of  Naval  Operations  (Op-32) 

Navy  Department 
Washington  25,  D.  C. 

Director 

Naval  Ordnauce  Laboratory- 
White  Oak,  Maryland 

Commande  r 

U.  S.  Naval  Electronics  Laboratory 
San  Diego,  California 

Commander  (AAEL) 

Naval  Air  ueveiopment  Center 
Jchnsville,  Pennsylvania 

Libra  rian 

U.  S.  Naval  Post  Graduate  School 
Monterey,  California 

Director 

Signal  Corps  Engineering  Laboratories 
Evans  Signal  Laboratory 
Supply  Receiving  Section 
Building  No.  42 
Belmar,  New  Jersey 


t 

w 


r 


3 Commanding  General  (RDRRP) 

Air  Research  and  Development  Command 
Post  Office  Box  1395 
Baltimore  3,  Maryland 

2 Commanding  General  (RDDDE) 

Air  Research  and  Development  Command 
Post  Office  Box  1395 
Baltimore  3,  Maryland 

1 Commanding  General  (WCRR) 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,  Ohio 

1 Commanding  General  (WCRRH) 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,  Ohio 

1 Commanding  General  (WCRE) 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,  Ohio 

2 Commanding  G c nc  ral(WCRET) 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,  Ohio 

1 Commanding  General  (WCREO) 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,  Ohio 

2 Commanding  General  (WCL.R) 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,  Ohio 

1 Commanding  General  (WCLRR) 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,  Ohio 

2 Technical  Library 
Commanding  General 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,  Ohio 

1 Commanding  General  (R.CREC-4C) 

Rome  Air  Development  Center 
Griffiss  Air  Force  Base 
Rome,  New  York 

1 Commanding  General  (RCR) 

Rome  Air  De.  .'lopment  Center 
Griffiss  An  Force  Base 
Rome,  New  York 

-iii  - 


r 


n 


2 Commanding  General  (RCRW) 

Rome  Air  Development.  Center 
Griffiss  Air  Force  Base 
Rome , New  York 

6 Commanding  General  (CRR) 

Air  Force  Cambridge  Research  Center 
230  Albany  Street 
Cambridge  39,  Massachusetts 

1 Commanding  General 

Technical  Library 

Air  Force  Cambridge  Research  Center 
230  Albany  Street 
Cambridge  39,  Massachusetts 

2 Director 

Air  University  Library 

Maxwell  Air  Force  Base,  Alabama 

1 Commander 

Patrick  Air  Force  Base; 

Cocoa,  Florida 

2 Chief,  Western  Division 

Air  Research  and  Development  Command 
P.  O.  Box  2035 
Pasadena,  California 

1 Chief,  European  Office 

Air  Research  and  Development  Command 
Shell  Building 
60  Rue  Ravenstein 
Brussels,  Belgium 

1 U.  S.  Coast  Guard  (EEE) 

1300  E Street,  N.  W. 

Washington,  D.  C. 

1 Assistant  Secretary  of  Defe  se 

(Research  and  Development; 

Research  and  Development  Board 
Department  of  Defense 
Washington  25,  D.  C. 

5 Armed  Services  Technical  Information  Agency 

Document  Service  Center 
Knott  Building 
Dayton  2,  Ohio 


-iv- 


L 


Director 

Division  14,  Librarian 

National  Bureau  of  Standards 

Connecticut  Avenue  and  Van  Ness  St.  , N.  W. 

Director 

Division  14,  Librarian 

National  Bureau  of  Standards 

Connecticut  Avenue  and  Van  Ness  St.  , N.  V/. 

Office  of  Technical  Services 
Department  of  Commerce 
Washington  25,  D.  C. 

Commanding  Officer  and  Director 
U.  S.  Underwater  Sound  Laboratory 
New  London,  Connecticut 

Federal  Telecommunications  Laboratories,  Inc  . 

Technical  Library 

500  Washington  Avenue 

Nutley,  New  Jersey 

Librarian 

Radio  Corporation  of  America 
RCA  Laboratories 
Princeton,  New  Jersey 

Sperry  Gyroscope  Company 
Engineering  Librarian 
Great  Neck,  L.  I.,  New  York 

Watson  Laboratories 
Library 

Red  Bank,  New  Jersey 

Professor  E.  Weber 
Polytechnic  Institute  of  Brooklyn 
99  Livingston  Street 
Brooklyn  2,  New  York 

University  of  California 
Department  of  Electrical  Engineering 
Berkeley,  California 

Dr.  E.  T.  Booth 
Hudson  Laboratories 
145  Palisade  Street 
Dobbs  Ferry,  New  York 

Cornell  University 

Department  of  Electrical  Engineering 
Ithaca,  New  York 


- v- 


Unive*  sity  of  Illinois 

Department  of  Electrical  Engineering 

Urbane. , Illinois 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
Silver  Spring,  Maryland 

Professor  A.  von  Hippel 
Massachusetts  Institute  of  Technology 
Research  Laboratory  for  Insulation  Research 
Cambridge,  Massachusetts 

Director 

Lincoln  Laboratory 

Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts 

Signal  Corps  Liaison  Office 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts 

Mr.  Hewitt 

Massachusetts  Institute  of  Technology 
Document  Room 

R.e search  Laboratory  of  Electronics 
Cambridge,  Massachusetts 

Stanford  University 
Electronics  Research  Laboratory 
Stanford,  California 

Professor  A.  W.  Straiton 
University  of  Texas 

Department  of  Electrical  Engineering 
Austin  12,  Texas 

Y ale  .University 

Department  of  Electrical  Engineering 
New  Haven,  Connecticut 

Mr.  James  F.  Trosch,  Administrative  Aide 
Columbia  Radiation  Laboratory 
Columbia  University 
538  West  120th  Street 
New  York.  27,  N.  Y . 

Dr.  J.V.N.  Granger 
Stanford  Research  Institute 
Stanford,  California 


-vi- 


